Quantum kinetic description of Coulomb effects 
in one-dimensional nano-transistors 
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In this article, we combine the modified electrostatics of a one- dimensional transistor structure 
with a quantum kinetic formulation of Coulomb interaction and nonequilibrium transport. A multi- 
configurational self-consistent Green's function approach is presented, accounting for fluctuating 
electron numbers. On this basis we provide a theory for the simulation of electronic transport and 
quantum charging effects in nano-transistors, such as gated carbon nanotube and whisker devices and 
one-dimensional CMOS transistors. Single-electron charging effects arise naturally as a consequence 
of the Coulomb repulsion within the channel. 



I. INTRODUCTION 

As scaling of field-effect-transistor (FET) devices 
reaches the deca-nanometer regime, multi-gate architec- 
tures and ultrathin active channel regions are manda- 
tory in order to preserve electrostatic integrity. It has 
been shown that a coaxially gated nanowire represents 
the ideal device structure for ultimately scaled FETsi^ 
A variety of ID nanostructures - such as carbon nan- 
otubes, silicon nanowires or compound semiconductor 
nano-whiskers - have been demonstrated and intensive re- 
search has been devoted to the realization of field-effect- 
transistor action in these nanostructuresii^*& Due to the 
small lateral extent in the nanometer range, electronic 
transport through such nanowires is one-dimensional 
with only a few or even a single transverse mode par- 
ticipating in the current. As a result, increasingly less 
electrons are involved in the switching of a nanowire tran- 
sistor. In fact, even in devices with rather long channel 
lengths of lOOnm, only on the order of 1-10 electrons 
constitute the charge in the channel for on-state volt- 
age conditions. Hence, single-electron charging effects 
are increasingly important and have to be taken into 
accounti 7 i 8 i 9 

Two different approaches are commonly used to de- 
scribe electronic transport in nano-transistors: A quan- 
tum kinetic approach based on real-time Green's func- 
tions provides an excellent description of non-equilibrium 
states^*^*^ Here, the Coulomb interaction is described 
in terms of a selfconsistent Hartree potential, option- 
ally combined with a spin-density-functional exchange- 
correlation term in local density approximation (LDA- 
SDFT). However, this framework does not account for 
single-electron charging effects without forcing integer 
electron numbers. Alternatively, the second approach 
considers a quasi-isolated nanosystem with a many-body 
formulation of Coulomb interaction, including electronic 
transport on a basis of rate-equationsii&iiiSiiSiiLiSiiii 
While predicting single-electron charging effects cor- 



rectly, the latter neglects dissipation and renormalization 
effects due to the source and drain contacts. 

Here, we present a novel approach that allows to com- 
bine a quantum kinetic description of non-equilibrium 
electron transport with non-local many-body Coulomb 
effects in one-dimensional FET nanodevices. Within our 
approach, single-electron charging effects arise naturally 
as a consequence of the Coulomb interaction. Our for- 
malism contains two central ingredients: In order to 
cope with particle-number fluctuations under nonequi- 
librium conditions, we introduce a multi-configurational 
self-consistent Green's function algorithm. Secondly, we 
consider a one-dimensional Coulomb Green's function for 
the transistor channel that allows to properly incorporate 
many-body interaction effects into a quantum kinetic ap- 
proach with electrostatic boundary conditions for a real- 
istic FET. As an example, we calculate the transfer char- 
acteristics of a nanowire transistor with Schottky-barriers 
(SB) at the contact-channel interfaces. 



II. COULOMB GREEN'S FUNCTION 




source/drain 
contact 



FIG. 1: Schematic view of a ID FET geometry. (d ox and d c ;j 
denote the gate insulator thickness and the channel diameter of 
the nano-transistor, respectively.) 



Consider a coaxially gated field-effect-transistor as il- 
lustrated in Fig. ^ A cylindrical semiconducting channel 
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material is surrounded by a thin dielectric and a metal- 
lic gate electrode. The electrostatic potential V inside 
such a one-dimensional (ID) transistor channel obeys a 
modified Poisson equation^^ 



dx 2 
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where p is the ID charge density. Vq denotes the gate 
potential and S = ird 2 ch /A is the effective cross-sectional 
area. The characteristic length A is related to the spa- 
tial separation of the gate electrode from the channel 
(which should be smaller than the total length L of the 
channel) iSS Note that Eq. is an appropriate descrip- 
tion for coaxial as well as planar transistor geometries, 
differing only in the characteristic length A. In the fol- 
lowing, we assume perfect metallic source and drain con- 
tacts at the boundaries, yielding fixed-potential bound- 
ary conditions due to given chemical potentials within 
these contacts. 

A key ingredient in our formalism is the usage of a 
Coulomb Green's function for the description of charge 
interaction within the channel. This allows us to formu- 
late classical electrostatics (with boundary conditions) 
and many-body interaction between electrons on equal 
footing. The corresponding Coulomb Green's function 
of the gated channel (with < x, x' < L and vanishing 
potential at the boundaries 0, L) can readily be obtained 
as 



v(x, x') 
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(In contrast, if we considered open boundary conditions, 
we would obtain v{x,x r ) = (A/2) exp (— \x — x'\/X) in- 
stead.) For a given charge density p inside the channel, 
the potential thus reads 



V(x) 



1 



dx v(x, x) p(x) + V ext {x), (3) 



with the external potential contribution 



V ext {x) 



sinhf-^r^) sinh (%■) 

. \ X J v s + —rmV D (4) 
sinh (f ) sinh (i) 
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+ A^ 



jdx' v(x,x') V G , 



where Vg and Vd denote the contact potentials. 



III. SYSTEM HAMILTONIAN 

In this article, we make use of a tight-binding descrip- 
tion of the system, represented by a localized ID single- 
particle basis {4>j(x, a)} (where the single-particle index 



j represents a combined orbital, site, and spin multi- 
index.) The total system Hamiltonian H = Hq + H ee + 
Hs + Hd can be split into four parts. Hq contains all 
single-particle terms of the transistor channel: 



Ck, 



(5) 



h 3k = ~e dx <f>* (x, a) <j> k (x, a) V dop (x) + V ext (x) 
+Sjkdj + tjk, 



with the electron annihilation operator Cj for state j. The 
composition of the channel (material-specific properties, 
layer sequence, etc.) is described by dj and off-diagonal 
hopping matrix elements t^iSi^ Vd op denotes the po- 
tential resulting from fixed charges pdop (due to ionized 
doping atoms), whereas V ex t stems from external charges 
due to the applied drain-source voltage and the gate in- 
fluence (see Eq. iJU). 

Furthermore, Hs and Hjj are the Hamiltonians for 
the source and drain contacts, respectively. Latter also 
contain the corresponding hopping terms to the outer 
ends of the channel region, providing electron injection 
and absorption. Each contact is assumed to be in a state 
of local equilibrium with an individual chemical potential 
according to the applied voltage. (See also Eqs. ©, i(T2"|) 
below.) 

Most importantly, H ee describes the many-body 
Coulomb interaction between electrons within the chan- 
nel region*^ 



with Coulomb matrix elements 

^2 



V, 



mjkl 



— -> Idxldx v(x.x') 
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(7) 



x <Kn( x i°') ^K^' ") <t>k(x,<x) <j>l(x,a), 

which employ the Coulomb Green's function Eq. @. 

IV. QUANTUM KINETICS 

A quantum kinetic description of the system (un- 
der nonequilibrium conditions in particular) is ob- 
tained via the usage of a real-time Green's functions 
formalismSSi^ The retarded and lesser (two-point) 
Green's functions in the time domain are given by 



cr jk (t) = -ie(t)<{ C j(t), 4(o)}>, 



(8) 



*<4(0)c*(t)>, 



for steady-state conditions. In the following, we will 
consider the Fourier transformed functions, defined via 
G{E) = (l/h)fdtexp(iEt/h)G(t). 
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For temperatures T well above the Kondo temper- 
ature of the system, the Coulomb interaction can be 
treated independently of the contact coupling, albeit self- 
consistently. In matrix notation, the Dyson equation for 
the channel thus can be written asi2i2£i2& 

G r = G r0 + G r0 [Y, r ee + Y> r s + E r D } G r , (9) 
G< = i f A 

+i (fs - fo) G r T s G a 

+i (Id - fo) G r T D G a , 

with T s = i(E r s - T D = i(E r D - H a D ) and A = 
i(G r — G a ). fs and fo are the local source and drain 
Fermi distribution functions, respectively, assuming lo- 
cal equilibrium within these reservoirs. On the other 
hand, fo denotes the equilibrium distribution function 
of the isolated channel system (typically, fo = fr>). Fur- 
thermore, G r0 = {E - h + ie)- 1 (with e -> 0+) and 
G a = G rt . 

Once G < has been obtained selfconsistently from 
Eq. 10, observables like the electron density p e and the 
current I e (through an arbitrary layer at xo) can be cal- 
culated via 

p e (x,a) = -e^2 <j>j(x,a) <f> k (x,<r) fa, (10) 
jh 

I e = -- ^2 2 Im(t jk pjk), 

xj<xo, 

with the single-particle density-matrix 

Pj k = ±-JdE\ G <{E). (11) 

The effective contact selfenergies due to the coupling 
of the channel to the source and drain regions (c = S, D) 
can be obtained a q 10 i 26 i 27 

^cjk( E ) = X! tc JP G c° P q( E ) tcqk, (12) 

with the isolated contact Green's function G£° and 
contact-channel hopping terms t c . 

The evaluation of the Coulomb selfenergy Y7 ee re- 
quires a suitable approximation scheme due to the in- 
finite Green's function hierarchy (which is a consequence 
of the two-particle interaction) . As a first-order expan- 
sion (Hartree-Fock diagrams), four-point Green's func- 
tions can be factorized into linear combinations of prod- 
ucts of two-point functions Using this approxima- 
tion, the Coulomb selfenergy reads 

S Lmi = X! (VmM ~ Vjmkl) Pjk- (13) 

Note that T,f e = 0, and Y7 ee is non-local, hermitian and 
energy-independent (static) within the considered ap- 
proximation scheme; compare also with Ref.l26i For con- 
venience, the Hartree potential (first V term in Eq. (113(1 ) 

V„(x) = — !— V fdx' v(x,x') Pe (x',a') (14) 

tO^chd , J 

rr' 



can be separated from the retarded Coulomb selfenergy 
(compare Eq. ©), where the electron charge density p e is 
given by Eq. (JTOJl. Hence, the total electrostatic potential 
of the system reads V = Vd op + Vh + V ex t ■ 

For integer-number electron filling conditions, Eq. 1|13|) 
provides an excellent description of the system for 
application-relevant temperatures. However, under 
nonequilibrium conditions, one has to deal with non- 
integer average filling situations, which are beyond the 
scope of a first-order (mean-field) selfenergy in general. 
In the following section, we will therefore present a multi- 
configurational approach which is able to cope with such 
particle-number fluctuations. 

V. MULTI-CONFIGURATIONAL 
SELF-CONSISTENT GREEN'S FUNCTION 

A thermodynamic state of the transistor channel with 
fluctuating electron number can be considered as a 
weighted mixture of many-body states with integer fill- 
ing (configurations) of relevant single-particle states. For 
a given G < , relevant single-particle states are defined as 
eigenstates of p (Eq. dU) that exhibit significant occu- 
pation fluctuations and have a sufficiently small dephas- 
ing (due to the contact-coupling). This projection to a 
relevant single-particle subspace of dimension N reduces 
the resulting Fock subspace dimension 2 N significantly, 
rendering this approach numerically feasible. 

For each configuration, the Coulomb selfenergy ap- 
proximation Eq. (|T3")l becomes adequate. Then the 
Green's function can be written as a configuration- 
average: 

G = ^> K G[p K ], (15) 

K 

where p K denotes the single-particle density-matrix (de- 
rived from p) for configuration k with weight w K . G[p K ] is 
the corresponding Green's function (retarded and lesser) 
which is obtained by using Eqs. ©, l|13|l . 

The weight vector w defines a projected nonequilib- 
rium many-body statistical operator in the relevant Fock 
basis. (Note that the configurations defined above might 
not be exact eigenstates of the projected many-body 
Hamiltonian, containing Eq. (jSJ in particular. In the 
following, we restrict ourselves to the dominant diago- 
nal elements of the many-body Hamiltonian in the rel- 
evant Fock basis.) Conscquentcly, the resulting many- 
body lesser Green's function reads 

GusHkit) = * 5>« e*<*-*»>* (16) 

K, A 

x (k|cJ,|A) (A|cj|k;), 

where \n) denotes a relevant Fock state with en- 
ergy E K . In principle, w must be chosen such that 
A(G^ /B [w],G < ) = min for a given G < within the rel- 
evant subspace, where A measures the cumulative dif- 
ference of spectral weights of corresponding resonances. 
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However, for most applications it is sufficient to con- 
sider a vector uu that maximizes the entropy at an ef- 
fective temperature T*(T,Vg,Vd) under the (weaker) 
subsidiary condition that pmb[w] = p within the rel- 
evant subspace, where pmb[w] denotes the many-body 
result ^2iSSi22i (Under moderate bias conditions, it is jus- 
tified to assume T* w T.) 

In turn, G < from Eq. H15fl can be taken as a new G < , 
serving as an input for the calculation scheme described 
above. This defines a self-consistency procedure for G 
and w, which we refer to as the multi-configurational 
self-consistent Green's function algorithm (MCSCG). 
Such an approach resembles the multi-configurational 
self-consistent field (MCSCF) approximation. 30 How- 
ever, MCSCG deals with grand-canonical noncquilibrium 
states and considers an incoherent superposition (mix- 
ture) of different configurations. Obviously, coherent 
superpositions of many-body states of varying particle 
numbers would be subject to strong dephasing due to 
the Coulomb interaction and the resulting entanglement 
with the environment. 

Having solved the many-body diagonalization problem 
of relevant states, it is straight-forward to employ this 
approach to calculate higher-order correlation functions 
(within the relevant subspace). Note that the MCSCG 
can also be interpreted as a means to construct a non- 
static S ee for Eq. 



VI. EXAMPLE: SB-FET 

In the following illustrative example, we consider a 
one-band nanowire-FET with Schottky-barrier injectors, 
having one localized orbital (with two spin directions) 
per site. Therefore, only Coulomb matrix elements of 
the form Vijji are remaining. Furthermore, we assume 
nearest-neighbor hopping with a real hopping parame- 
ter t = Ti 2 / {2m* a 2 ) . We have used the following device 
parameters: The nominally undoped channel has a di- 
ameter of d c h — 4nm and a length of L = 20nm (im- 
plemented as 20 sites with a spacing of a = lnm). The 
channel with e c h = 15 is surrounded by a gate oxide with 
dox — lOnm and e ox = 3.9, yielding A ~ 3.7nm. We as- 
sume an effective electron mass of m* = 0.05m e (giving 
t = 0.77eV). The metallic source and drain contacts have 
a Schottky-barrier height of $ss = 0.5eV. For simplicity, 
the corresponding contact selfenergy is assumed to be of 
the form S c « — iF/2 (within the band at the outer ends 
of the channel) with T w 76meV. Note that this param- 
eter has to be chosen to match the actual metal contact 
used in an experiment. However, it is uncritical for the 
electronic spectrum and single-electron charging effects. 
The system temperature is T = 77K. Up to N — 6 adap- 
tively chosen relevant single-particle states were taken 
into account, depending on the applied voltages (with 
V s = 0). 

Fig. |21 visualizes the local density of states (LDOS) for 
low drain-source bias conditions and two different gate 
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FIG. 2: Spectral function A(x, E) as a grayscale plot for (a) 
V G = 0.59V and (b) V G = 0.71V. In both cases, V D = 2mV is 
chosen. The resulting average electron number within the channel 
is (a) N e ft: and (b) N e rj 1 (the electron is located in the 
lowest resonance). The solid white line represents the mean-field 
potential V(x), whereas fig and /ip denote the chemical potentials 
of the source and drain contact, respectively. T = 77K. 



voltages Vq = 0.59V and Vq = 0.71V, where the av- 
erage electron number in the channel becomes N e 
and N e 1, respectively. The existence of quasi-bound 
states (i.e. spatially and energetically localized reso- 
nances in the spectral function A) yields discrete single- 
electron energies with associated interaction energies due 
to Sg e . Comparing the situation for N e = 1 with jV e = 0, 
the single-electron resonances are moved to higher ener- 
gies with respect to the lowest energy state due to the 
Coulomb repulsion. Note that each electron is not sub- 
ject to its own Hartree potential (see lowest resonance 
in Fig. [2Kb)) because Y7 ee does not contain unphysical 
self-interaction energies, but includes exchange terms and 
correctly accounts for the electron spin. In the shown ex- 
ample, the next higher available state for a second elec- 
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tron (with opposite spin) is separated by the interaction 
energy Voo ~ 93meV (see arrow in Fig. Hfb)). In gen- 
eral, energy levels are splitted by exchange energy terms, 
which have a significant influence on the energy spec- 
trum. 

As a natural consequence of ft. + Y7 ee we therefore ex- 
pect to observe the effect of a step-like electron filling 
(under conditions close to equilibrium in particular), en- 
ergetically determined by single-electron levels and re- 
pulsion energies. This behavior in fact can be seen in 
Fig. |3 where the electron filling characteristics is plot- 
ted for a varying gate voltage and fixed drain-source bias 
Vd = 2mV. Furthermore, Coulomb oscillations in the 
accompanying current through the channel can be iden- 
tified. 



j i i i i i i i i i i i i i i i i_ 
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GATE VOLTAGE (V) 

FIG. 3: Single-electron tunneling characteristics for Vd = 2mV. 
The solid line with filled squares corresponds to the average elec- 
tron number within the potential well, whereas the dashed line with 
open circles shows the drain current, exhibiting Coulomb oscilla- 
tions. T = 77K. 

Models solely based on a selfconsistent Hartree po- 
tential do not provide such quantization effects due to 
Coulomb repulsion. With a Hartree approach (as used 
with conventional Schrodinger-Poisson solvers), spectral 
features are solely shifted in energy, depending on the av- 
erage electron density. In contrast, the MCSCG (as well 
as the exact diagonalization of the isolated system) pro- 
vides a superposition of fading spectra of integer electron 
numbers with full interaction energies, however, having 
spectral weights that depend on the average filling con- 
dition. The local density of states under nonequilibrium 
conditions as shown in Fig. 0] clearly demonstrates this 
behavior, where the average electron number within the 
well is iV e w 0.22. In fact, the expectation value of the 
electron number need not be an integer, especially under 
non-equilibrium bias conditions, which can be seen in 
the corresponding transfer characteristics of the system 
as plotted in Fig. 03 Furthermore, Fig. visualizes the 
output IV characteristics, where a finite drain voltage is 
required to pull the chemical potential of the drain con- 
tact below the lowest energy level. These results clearly 
demonstrate the strengths of the MCSCG approach, be- 
ing able to describe single-electron charging effects under 
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FIG. 4: Nonequilibrium spectral function A(x, E) as a grayscale 
plot for Vq = 0.7V and Vd = 0.2V. The resulting average elec- 
tron number within the channel is N e ~ 0.22. (The solid white 
line represents the mean-field potential V(x), whereas fig and (id 
denote the chemical potentials of the source and drain contact, 
respectively.) T = 77K. 
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FIG. 5: Transfer characteristics for Vd = 0.2V. The solid line with 
filled squares corresponds to the average electron number within the 
potential well, whereas the dashed line with open circles shows the 
drain current through the channel. T = 77K. 



nonequilibrium bias conditions with fluctuating electron 
numbers. 

In general, we expect the many-body Coulomb interac- 
tion to have a significant impact on the electrical behavior 
of nano-transistors if the single-electron charging energy 
becomes > 4fcT, having consequences for the transcon- 
ductance, onset/pinch-off voltages, sub-threshold cur- 
rents, and system capacitance. A more detailed discus- 
sion of these aspects will be published elsewhere. 
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FIG. 6: Output characteristics for V G = 0.59V and V G = 0.71V. 
T = 77K. 



of electronic transport enables us to describe many-body 
charging effects within the transistor channel. We have 
presented a multi-configurational self-consistent Green's 
function algorithm, which is able to cope with fluc- 
tuating electron numbers under noncquilibrium condi- 
tions. In the discussed example of a nano-FET with 
Schottky-barrier injectors, we have visualized how single- 
electron charging effects arise naturally as a consequence 
of the many-body Coulomb repulsion between quasi- 
bound states. The usage of a Green's function formula- 
tion permits the systematic extension to further Coulomb 
diagrams and the consistent inclusion of phonon scatter- 
ing. 



VII. CONCLUSION 

The Coulomb Green's function of a one-dimensional 
FET in combination with a quantum kinetic description 



With the presented theoretical approach, we are able 
to describe electronic transport and quantum charging 
effects in ID nano-transistors such as gated carbon nan- 
otubes, semiconductor whiskers, and ID CMOS transis- 
tors (in coaxial and multi-gate planar geometry). 
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